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ABSTRACT 


Differential corrections analyses of three W UMa type binaries show that 
RZ Tau and AW UMa have common envelopes and the relatively small gravity 
darkening predicted by Lucy, while R Z Comae seems to have a larger gravity 
effect and is only marginally in contact. Some details of our method for com- 

f 

puting contact binary light curves are given. We suggest that W UMa binaries 
which have occultati on primary eclipses may, in general, have large amplitudes 
for gravity darkening or a photometric surface brightness effect which simulates 
large gravity darkening. 



INTRODUCTION 


The term contact binary has been in use for more than two decades, 
but until recently there has been some uncertainty as to how literally 
this designation should be taken. That is, the traditional methods for 
light curve analysis, which are based on the Russell model of similar 
ellipsoids, break down for very close binaries and cannot be used to 
compute the light variation for components which are fully in contact. 
Recently Lucy (1967, 1968a, 1968b), Rucinski (1969), Whelan (1972), 
Biermann and Thomas (1972), and others have put forward interesting 
ideas on the nature of contact binaries which are subject to quantita- 
tive test against the available observations, provided that suitable means 
are developed to compute the light curves. Such tests have been made 
by Lucy (1971) and by Mochnacki and Doughty ' (l 972a, 1972b) who found that: 

a) At least some. 1 of the W UMa type binaries are fully in contact 
(i.e. they have common envelopes whose "photospheres" follow a single 
equipotential surface). 

b) The light variation for some contact binaries is compatible with 
the gravity darkening law predicted by Lucy for convective envelopes, 
which has an amplitude about one-third that of Von Zeipel (1924) radiative 


gravity darkening. 
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Sixteen systems were analyzed by Lucy and six by Mochnacki and Doughty. 
In this paper we give results of our analyses of only three systems (RZ 
Comae, RZ Tauri, and AW UMa) but these have been fitted by a procedure 
which has significant advantages over those by Lucy and by Mochnacki and 
•Doughty. Inthe major points, our conclusions are similar to those by 
the other authors, but we have found the gravity darkening law by least 
squares, while they have assumed the "convective' 1 law to hold. We list 
below the reasons we consider quantitative results based on our approach 
(Wilson and Devinney, 1971) to be relatively well determined. 

a) The accuracy of the basic light curve program is at least one 

order of magnitude better than that of Lucy or Mochnacki and Doughty, 
who claim accuracy of only about ±0 .01. This is important for dif- 

ferential corrections solutions since the required numerical derivatives 
are computed by differencing adjacent values of the dependent variable. 

The light curves differenced in our solutions have a precision of one 

part in several thousand in the worst case or, more typically, one part in 

io 4 . 

b) The results, whibh we give in Table 1, follow from least 
squares solutions which make use of observations in two passbands 



{for HZ Tau and AW UMa, but not for RZ Com) fitted simultaneously. 

In fact, an arbitrary number of passbands can be treated simultaneously 
but we presently avoid ultraviolet light curves because of blanketing 
effects. Observation points are weighted by the program according to the 
source of noise (variable sky transparency, shot noise, instrumental noise); 
according to the relative scatter of the several light curves; and according 

to the intrinsic weight of individual points (eg. number of individual obser- 

*\ 

vations per normal point). For 3-11 three binaries discussed here we 

■. # i 

assumed variable transparency to be the main source of scatter. 

c) The fit to the observations seems quite good for all three binaries, 
with some minor qualifications, which will be noted. 

d) A very stringent test of the degree to which the model represents 

a real binary is afforded by including third light as an adjustable parameter. 
Models which are physically unrealistic can be fitted in many cases by 
•'trading- off" their deficiencies for third light. When third light was included 
as a free parameter, we found it to be of the order of its probable error. 

e) Before applying the differential corrections program, reasonably 
extensive trial-and-error experiments were carried out with the light curve 
program. For each binary, therefore, even the first iteration of the differ- 
ential corrections program began from a relatively good fit. 

Incidentally, the differential corrections program is now faster than that 


used for MR Cyg (Wilson and Devinney, 1971) by a factor of 8. 



MODEL CONSTRAINTS FOR CONTACT BINARIES 


The basic model and computing scheme used here has been des- 
cribed by Wilson and Devinney (1971), and some additional details of the 
reflection effect have been given by Wilson, et. al. (1972), so we need 
only state some constraints imposed when operating in the contact binary 
mode. Since our program treats the two binary components separately, we 
specify that there be no physical discontinuities across the boundary (on 
the common envelope) which separates, for computing purposes, one com- 
ponent from the other. To some extent the placement of this boundary is 
arbitrary, and we define it to be the intersection of the common envelope 
with a plane which passes through the inner Lagrangian point normal to 
the line of centers, or x-axis. This boundary curve is therefore a ring 
around the narrow part of the neck connecting the components, and in- 
tegration of the light of a component begins and ends at this ring for 

i 

the lower latitude rows on each component. The parameter list for con- 
tact binaries includes: 

i : inclination of the orbit plane to the plane of the sky, 
g : exponent in the gravity darkening law. 



: polar effective temperature for the component covered at 
primary eclipse (component l) 

A : bolometric albedo, 

Q : Roche modified surface potential, 

q : mass ratio, m /m . 

£ I 

1^ : relative luminosity of component l, 

x : limb darkening coefficient 

^2 : third light t 

Additional details on our conventions regarding g, Q, q, L^, and x were 

given by Wilson and Devinney (1971). Complete definitions for A and t 
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are to be found in Wilson, et. al. (1972). When these parameters have 
been specified, the polar temperature of the secondary is constrained to 
have that value which causes the boundary discontinuity in local surface 
flux to vanish. This temperature may be computed from the gravity 
darkening law as 


T 2 = T I (VO 


0.25 g 
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where a^ a^ are quantities proportional to the acceleration due to gravity 
at the poles of the two components (this notation is unconventional, but we 
have been using g for the gravity darkening exponent). The luminosity 
of the secondary component, L 2 » is now determined. L 2 is computed by 



numerical integration but is not, of course, a free parameter, so we 
list only L (normalized) = L / (L + L ). Indeed, L is a free parameter 

i i i £ L 

only in the sense of being a convenient scale factor (before normalization), 
since the luminosit}' ratio is entirely determined for our contact binary 
model when i, g, T^, 0, and q have been fixed. To complete the list of 

= 0^, g 2 = A 2 = A Jf and x 0 = x^ automatically 
when running in the contact mode. It may be desirable in the future to 
replace the last three conditions by more sophisticated constraints based 
on improved physical models for contact binaries and their energy exchange. 

Our approach differs from that of Mochnacki and Doughty and of Lucy 
in that we treat the limb and gravity darkenings as adjustable parameters. 
It would seem that both approaches can be defended. We feel that, at 
this early stage, departures of these parameters from expected values 
may offer some of the best clues regarding necessary modifications for 
theoretical models. Also, in view of the demonstrated strong determinacy 
of photometric mass ratios, Lucy’s fixing of q at the (sometimes decidedly 
incorrect) spectroscopic value might vitiate some solutions to the extent 
that the deduced degree of contact might be seriously affected. 


constraints, we set 
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RESULTS 

The systems analyzed here were selected on the basis of having 
accurate observations, complete eclipses, and absence of light curve 
assymetries. AW UMa does, in fact, have a light curve which is assymetric 
about phase zero, but it is extremely interesting from certain points of 
view and was therefore included. All three binaries were analyzed by Lucy 
(1971), while RZ Tau and AW UMa were studied also by Mochnacki and 
Doughty (1972a, 1972b). The final adjusted parameters from our differen- 
tial corrections program are listed in Table l and the resulting theoretical 
curves are graphed among the observations in Figures 1, 2, and 3. Also 
given in Table l are absolute dimensions and masses found by incorpor- 
ating radial velocity information (Paczynski, 1964; Struve and Gratton 1948). 
Absolute dimensions are not listed for RZ Tau because the spectroscopic 
observations (Struve, et. al. , 1950) are not sufficiently numerous to yield 
meaningful results. In computing these quantities we used our photometric 
mass ratios in combination with spectroscopic information from the 
radial velocity curve of the more massive component. 

For both RZ Com and AW UMa we find masses which seem incom- 
patible with their positions in the HR diagram. The more massive component 
of each binary is only slightly above the main sequence but has a mass appro- 
priate to a much earlier star. RZ Com, for example, seems to have the 
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mass (2.0 m ) of an early A main sequence star, yet is only about one magnitude 
above the main sequence at type K0. If the spectroscopic rather than the photo- 
metric mass ratio is used, we find 1.6 m rather than 2.0 m , which is still 

G 0 

too high and really does not resolve the difficulty. Taken at face value, the 
masses for RZ Com and AW UMa are not consistent with any normal view of 
pre or post-main sequence evolution, nor even' with any reasonably simple 
binary star evolution. However, according to Popper (1967), radial velocity 
amplitudes are measured systematically too large from plates taken at inade- 

t 

quate dispersion (> 20 A/mm). The dispersions for RZ Com and AW UMa were, 
respectively, 76 and 75 A/mm at H y. Furthermore Paczynski (1964) speci- 
fically commented that he suspected blending problems for AW UMa, so the 
main problem seems to be that for RZ Com, whose radial velocity curves seem 
to be reasonably well defined. Probably this system should be observed again 
at high dispersion. 

For the binaries with observations in two colors we find, as usual, one 
value for each of the elements i, g, 0, and q; and two values (for the two 
passbands) for and x. RZ Comae is represented by two significantly 
different X5300 A light curves observed (Broglia, 1960) at slightly different 
epochs, which we analyzed separately. The resulting solution differences 

i 

for the two epochs will be discussed. 
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Table l also contains theoretical values for linear cosine limb darkening 
coefficients, determined by graphical fitting to tables by Carbon and 
Gingericli (1968), for comparison with our observational values. The. agree- 
ment seems satisfactoi'y considering that these numbers are probably quite 
sensitive to small photometric. perturbations and deficiencies of the model. 


The bolometric albedo. A, when allowed to adjust, came in most 
cases to values near unity, but with rather large probable error. Indeed, 

i 

one can expect the albedo to be weakly determined because the reflection 
effect is fairly small for contact systems. Although Rucinski (1969) has 
predicted that A should lie in the range 0.4 to 0. 5 for stars with fully adia- 
batic convective envelopes, and this has been verified by Hosokawa (1957, 1959), 
Napier (1970), and Wilson, et. al. ( 1972 ) for Algol and other binaries, we have 
simply set A to unity for all final iterations. We do not claim that is the 
ideal way to treat this element - only that it is one of several possible ways 
and the one we actually followed. Since the reflection effect is small, effects 
on the adjustment of other parameters should be slight. 

Mochnacki and Doughty (1972b) have pointed out an apparent incon- 
sistency between the spectral type and color index (Eggen, 1967) of RZ 
Tau, and suggested that the spectral type is no earlier than F6, rather 
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than FO as given by Struve, et. al. (1950). We have assumed that the spectral 
type is correct since it is not yet established that interstellar reddening is not 
responsible for the discrepancy. According to Eggen (1967), the color excess 
and distance modulus are nearly indeterminate for this system. In any case, 
the solutions are very insensitive to changes in the assumed temperature of the 
primary component. Effective temperatures for all three systems were estimated 
from the Harris (1963) calibration of spectral type against T . These tempera- 
tures were made slightly higher than the Harris values because they are polar ' 
temperatures in the program — not mean surface temperatures. 
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DISCUSSION 

AVV-UMa - The light curves of this system were observed by Paczynski 
(1964). It has a significant assymetry but there is no effect on the listed para- 
meters due to this problem. We repeated die solution after having removed 
the assymetry (by subtraction of a sine wave) and found essentially the same 
results. The final computed light curves "split the difference" between die light 
levels at the two quadrature points. There is a small discrepancy in secondary 
eclipse between the fits in B_ and V, but we do not consider diis alarming in 
view of the number of constraints enforced in the simultaneous two-color 
solutions, and consequent small number of degrees of freedom. Our values for 
the inclination, degree of contact and the rather extreme mass ratio are similar 
to those found by Lucy and by Mochnacki and Doughty, and the value for the 

_ o 

gravity darkening exponent, g = 4 j3 = 0.45 ± .06 provides impersonal evidence 
on the observational question posed by Lucy’s predicted value of 0. 32. The 
light curve assymetry rnigjht conceivably be caused by a stream of material 
which leaves the secondary component at its outer Lagrangian point if it could 
be shown that such a stream would have sufficient continuous opacity to atten- 
uate the light of the binary by about one percent. However we, along with Lucy 
and Mochnacki and Doughty, find the photosphere to be definitely within the 
outer contact surface, sb a positive ejection mechanism is required if this 


idea is to work. 
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RZ Tauri - Fitting of these observations (Binnendijk, 1963) gave little 
trouble after we abandoned the spectroscopic mass ratio (Struve, et. al. 1950) 
of 0.54. The agreement at X4420 and X'5300 seems very good except in the 
annular phases of primary eclipse, where we cannot reproduce the observed 
flat phase interval without destroying the fit at other phases. Mochnacki and 
Doughty had this same difficulty. We speculate that perhaps one or two un- 
fortunate residuals are responsible for this apparent flat region, and that it 

. A 

may not be shown by future light curves of RZ Tau. The value of g{0.29 ± 

. 06 p. e. ) is in excellent agreement with the convective value. A paper by one 
of us (Wilson, et. al. 1972) contains a remark that classical gravity darkening 
was used for the secondary component of the Algol system, based on experience , 
in fitting contact binaries. At that time we had completed our analysis of 
RZ Comae, for which we found g greater than unity, and were only beginning 
our study of the other two systems. These other systems, of course, did show 
the convective darkening. 

Insertion of the photometric mass ratio (0.37) into the mass function for 
the primary component gives individual masses of 4.0 and 1.5 which are 
almost surely too high. However one can see by plotting the few radial velocity 
points that the mass function itself is quite uncertain. A reduction in sin i 
of 25 percent would produce reasonable masses (about 1.6 and 0.6 with 
q = 0.37. It is probably a coincidence that the published values of sin i and 
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m /m happen to give plausible main sequence masses, since a small change 
2 1 

in either quantity, within the observational uncertainty, can give a result such 


as the above value of = 4. 0 m^. 

RZ Comae — The two X5300 A light curves shown in Figure 3 were ob- 
served two to three months apart by Broglia (1960). They are significantly 
different by inspection, and Broglia attributed this temporal behavior to a bright 
area existing on the outer facing hemisphere of the small star at epoch II 
(lower curve). Binnendijk, (1964) rediscussed these observations and pointed 
out that a somewhat more consistent explanation can be given in terms of a 
dark area on the inner hemisphere of the larger star at the other epoch. The 
decision between these possibilities is equivalent to a decision as to which 
epoch (if either) corresponds to the normal light curve of the system. Since 
the fit we have achieved is good at both epochs, one would tend to say that the 

light curve which yields parameters closer to the; expected (theoretical) values 

! 

is the normal case. The only elements for which we have predicted values are 
the limb and gravity darkening parameters. Unfortunately, the limb darkening 
is closer to the expected result for light curve II, while the gravity darkening is 
in better agreement for light curve I* Therefore it seems that we cannot, at 
present, say which of the two light curves is normal, and one might suspect 

t 

that transient effects are present at both epochs. 
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Unless these transient effects have seriously undermined the validity of 
our solutions for RZ Comae, it seems that the components are only marginally 
in contact. That is, one solution (epoch I) shows the binary to be slightly over 
contact while the other shows it slightly under contact. Could this result be 
related to the fact thatRZ Comae shows, apparently, a radiative (Von Zeipel) 
gravity effect, or perhaps an even larger effect? The spectral type of KO leads 
us to expect deep convection zones for the components of this system, whether 
they are in contact or not, and we do not know of any physical reason to associ- 
ate marginal contact with (apparently) radiative gravity darkening. However we 
can point to an interesting photometric effect, which may eventually lead us to 
understand the finding of a gravity exponent of unity for a convective envelope 
binary, when more systems with total primary eclipses have been studied 
sufficiently well. We have found, by experimenting with RZ Comae — like light 
curves, that the total eclipse can be the deeper only if g is of the order of unity 
or greater. No other parameter can be changed to produce this effect and, in 
particular, the convective value of g always yields annular primary eclipses. 

If these statements prove to be valid for all combinations of contact binary 
parameters, it will be easy to identify systems (they will be Binnendijk’s 
'class W) which, from the point of view of light curve analysis, show radiative 

i = 

gravity darkening. One can then hope to find common physical characteristics 
of these systems which mil indicate the reason for the observed darkening. 

For example, it may be found that such systems are subject to a particular 



IV 


photometric surface effect which simulates gravity darkening, although perhaps 
is physically imr elated to it. 

The last remarks serve to emphasize a point which many workers 
have recognized for some years - that since unexplained odd functions 
(sine terms) are present in the light curves of many systems, we should 
not be surprised to find even functions due to the same effects. An urgent 
present need, therefore, is for more quantitative work - both observational 
and theoretical - on the physical causes of assymetries in the light curves 
of close binaries. * ■ 
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Tabte i 


Parameters for Three Binaries 



RZ Coni (I) 

RZ Com tri) 

RZ Tau 

AW UMa 

i 

j 

SG° . 04 ± ?5i 

85° . 72 =?31 

82° .88 ± .35 

‘ 79°I0 ± .36 

e 

1. 13 x . 04 

l. 51 x . 02 

0.29 ± .06 

0.45 ± .06 

fcole) 

5500 K 

5500 K 

7200 K 

7000 K 

T 2 (pole) 

5564 K 

5552 K 

7146 K 

6737 K 

A 

LOO 

1. 00 

1.00 

1. 00 

n 

5. 618 ± . 054 

5. 869 ± .040 

2.4957 ± .0048 

i. 8321 ± . 0018 

q 

2.292 ±,030 

2.3 94 ± .020 

0.3721 x .0047 

0.0716 ± .', 005 

Lj &reen) 

0.3052 ± .0024 

0.2976 ±.0011 

0.7046 ± . 0057 

0. 9148 ±. 00 99 

L, (blue) 



0.7061 ±. 0062 

0. 9176 ± . 0125 

x (green) 

0. 50 ± ' . 03 

0. 62 ± . 02 

0.51 ±.03 

0.63 ±.01 

x (blue) 



. 57 ± . 03 

.69± .01 

Tj (pole) 

0.2 924 ± .0044 

0.2805 ± .0030 

0.4632 ±.0010 

0.5650 ± .OC06 

Tj (side) 

.3056 ±. 0052 

.2 918 ±.0035 

. 5024 ± . 0014 

. 6457 ± . 0011 

Tj (bank) 

.3403 ± .0082 

.3211 ± .0052 

5392 ±.0019 

.6644 ± . 0012 

r 2 feole) 

.4287 ±.0042 

.4240 ± .0029 

.3035 ±.00U 

.I860 ± . 0008 

r 2 (side) 

.4574 ± . 0055 

. 4509 ± .0037 

.3209 ±.0014 

.1956 ± .0010 

r 2 fcack) 

''.4859 ± .0071 

.4761 ± .0047 

.3802 ±. 0030 

.2537 ± ,0039 

Rj (side) 

0.63 X 10 6 km 

0. 60 X 10 6 km 


1.6 X 10 6 km 

R 2 (side) 

0. 95 X 10 km 

0. 93 X 10 6 km 


0.5X10 6 km 

m i 

0,9 m 

0 

0. 9 m 

0 


3. 1 m 

m 2 

2.0 m 

0 

2,1m 

0 


o 

0.2 m 

0 



Theoretical Values 



x (£reen) 

0.73 

0.73 

0. 66 

0, 65 

x (blue) 

•*•#»« 

i 

M 

0.77 

0.78 

0 (inner contact) 

5. 6599 

5. 8003 

2.6206 

l. 8670 

(outer contact) 

S. 0546 

5.1926 

2.3923 

1.8200 
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Captions for Figures 


Figure 1. 

Paczynski’s B, V observations of AW UMa and our theoretical curves. 
Figure 2. 

Binnendijk's U420 (blue) and >5300 green) observations of RZ Tauri with 
our fitted curves. Green i$ above. 

F igure 3 

Broglia's X5300 observations of RZ Comae at his epochs I (above) and II, 


with our theoretical curves. 
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